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1. Introduction 

Hard processes are ideal tools to characterize the medium produced in heavy-ion collisions []1|, 
At large enough virtuality the perturbative partonic cross section is unchanged and all medium 
effects appear as modifications of the long distance parton distributions (PDF) and fragmentation 
functions (FF). An excellent example is the production of particles at high transverse momentum 
in which large modifications of the fragmentation functions have been experimentally observed [Q. 

In the vacuum, the fragmentation of a perturbatively produced parton is well described by re- 
summing the leading soft and collinear singularities leading to an ordered shower structure described 
by the DGLAP evolution equations In heavy-ion collisions, an enhanced splitting probability 
is expected due to the additional medium-induced radiation by gluon bremsstrahlung off the fast 
parton, producing a softening of the associated jet structures. This effect is commonly identified as 
the dominant dynamical mechanism producing the strong suppression of high-p^ hadrons observed 
experimentally at RHIC 0. Provided the geometry of the collision is known, the only parame- 
ter controlling this suppression is the transport coefficient, q, characterizing the average squared 
transverse momentum acquired by the emitted gluon per mean free path in the medium. Hence, 
a one-to-one correspondence between the energy loss and the jet broadening exists as both are 
controlled by the same parameter 0. 

The above properties are based on the one-gluon inclusive spectrum computed in Refs. [|, |^, 0, |[ 
(see also related approaches in M, |TUi). For practical applications, however, more differential quan- 



tities, including exclusive one, two,. . . gluon emissions are needed. A successful phenomenology, 
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based on this radiative energy loss, assumes a separation between vacuum and medium contribu- 
tions to the shower development in which the "medium" radiation occurs before in time than the 
"vacuum" radiation, the latter taking place after the fast parton exits the medium. In this way, an 
independent emission approximation is taken for the first [|ll|, [12|, |13|, |14| while the usual (vacuum) 



DGLAP evolution describes the second. In the eikonal approximation for multiple gluon radiation, 
the energy loss distributions are Poissonian and normally known as quenching weights [0, |12|. One 
caveat of this formalism is the different role of the medium and vacuum radiations which looks 
artificial. Further limitations are the treatment of energy constrains and the role of virtuality. In 
the present paper we study these limitations and propose a new formalism in which the medium 
corrections are systematically included as an additive term in the splitting functions. This new 
term is directly taken from the known medium-induced gluon radiation spectrum by comparing 
the leading vacuum contributions. Our main assumption is the existence of an ordering variable 
also for the medium, which we will take to be the virtuality. This assumption is in agreement with 
the findings in and has been used before in |]T5|, |16|, |T^. Interestingly, we show that when the 
virtuality and the kinematic constrains coincide, the usual formalism, which provides an excellent 
agreement with experimental data on inclusive high-p^^ particle production, is recovered. 

Our more refined treatment of the fragmentation functions is mainly motivated by the immi- 
nence of the LHC heavy-ion program, in which real jets will be measured in heavy-ion collisions 



for the first time and their fragmentation functions reconstructed fl^, ^ . We present full 

calculations of medium-modified fragmentation functions at different energies and virtualities easily 
reachable within the first years of the program. Although our overall aim is the LHC physics, our 
work will also find applications to RHIC physics. In particular, we recalculate the high-p^- suppres- 
sion of light hadrons and check that the results agree with the ones already obtained within the 
Poisson approximation using the quenching weights. 

The paper is organized as follows: In Section |^ we first present the medium-modified splitting 
functions, which is then used to compute the Sudakov form factors and DGLAP evolution. In 
Section ^ we present our results for the medium- modified fragmentation functions and calculate 
the inclusive particle suppression to check, in Section ^ the degree of agreement with the previous 
formalism. Finally we present our conclusions. In Appendix ^ a derivation of the quenching weights 
is obtained as a limiting case of our formalism for large fraction of momentum and virtuality. 



2. Medium-evolved fragmentation functions 

The DGLAP evolution equations for the fragmentation function D{x, t) reads, omitting parton type 
indices. 

Here, P{z) is the splitting function describing the branching of a parton into two new ones with 
fractions of momenta z and 1 — z. For our implementation of medium effects in the fragmentation 
functions we recall the probabilistic interpretation of the (LO) DGLAP evolution which is evident 
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from its integral formulation (see e.g. p2| ) 



D{x, t) = A{t)D{x, to) + A(t) £ / ^ Pi^)D ti) . (2.2) 

The first term on the right-hand side in this expression corresponds to the contribution with no 
splittings between to and t while the second one gives the evolution when some finite amount of 
radiation is present. The evolution is controlled by the Sudakov form factors 



A{t) = exp 



to 



t' 27r ^ ' ^ 



(2.3) 



with the interpretation of the probability of no resolvable branching between the two scales t and 
to. 

The definition of the Sudakov form factors and its probabilistic interpretation depend on the 
cancellation of the different divergencies appearing in the corresponding Feynman diagrams, see e.g. 
p3| . Although such cancellation has never been proved on general grounds for partons re-scattering 
in a medium, it has been found in that, under certain assumptions, all the medium effects can 
be included in a redefinition of the splitting function 

P^°\z) = P-'^^iz) + AP{z, t), (2.4) 

where we have labeled as "vac" the corresponding vacuum splitting function. The main assumption 
to arrive at ( |2.4| ) is the independence of the multiple gluon emission when re-scattering with the 
nuclei is present so that an exponentiation of the splittings is possible. This independence could 
be broken if non-trivial color reconnections are present |2^, ^ but the size of this corrections for 
the DGLAP kinematics has not been computed. In this exploratory study, we will assume that 
the definition ( p. 4)) can be extended to the medium case. Similar assumptions have been done in 
15, 0, 1^ using simplified forms of the medium term AP which we will now improve by taking 



the full splitting probability as given by the medium-induced gluon radiation spectrum. As we will 
show, this procedure recovers the usually employed quenching weights [jTT], |12[ , making a connection 
with previous phenomenology in the field. 

We now describe the steps leading from the single inclusive distribution of gluons radiated 
in a medium from a parent parton, to the medium-modified fragmentation functions. We start 
by defining our medium-modified splitting functions, then we compute the corresponding Sudakov 
form factors, to finally write the modified DGLAP-like evolution equations that we use, and the 
initial conditions for such evolution. 

2.1 Medium-modified splitting functions 

The formalism of medium-induced gluon radiation p, |^, ^ provides the single inclusive distribu- 
tion of gluons emitted with energy u and transverse momentum k^, by medium-induced radiation 
from a parent parton with energy E, x = u/E. For the case of a massless parent the general formula 
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for the emitted energy reads 



dl 



a., C, 
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(2.5) 



Here, the Casimir factor Cr = | for a parent quark and 3 for a parent gluon, determines the 
couphng strength of the emitted gluon to the parent. Eq. ( p.5| ) re-sums the effects of arbitrary 
many medium-induced scatterings to leading order in 1/E. Properties of the medium enter (|2.5| ) 
via the product of the time-dependent density n(^) of scattering centers times the strength of a 
single elastic scattering o"(r). A detailed discussion of ( |2.5| ) including the physical interpretation of 
the integration variables {yi, yi, y, u, ^) can be found in Refs. 0, for the massless case and in 
p6| for the massive one. 

In ( |2.5|) the integration in yi [yi) signals the longitudinal position of the gluon radiation vertex 
in the amplitude (complex conjugate amplitude). For the case of high-p^ particle production, the 
parent parton is produced inside the medium and three cases can be distinguished depending on the 
position of the radiation vertex begin inside or outside the medium in both amplitude and complex 
conjugate amplitude and the interference of the two. In the case that both radiation vertex are 
outside the medium the usual vacuum radiation spectrum is obtained. This spectrum is normally 
subtracted from the total to define the medium-induced gluon radiation I^'^'^ as 



dl 



dp 



u- 



LU- 



+ UJ- 



dl 



med 



(2 6) 

du! dk± duj (ikx ' du (ik_|_ 
The vacuum spectrum defines the (vacuum) splitting function, which at small fraction x = 1 — z is 

(2.7) 



dp 



a 



dz dk^ 



1 



For the medium part, two different approximations are usually employed in ( p.5| ), the exact 
solution of the path integral being unknown. The first one corresponds to expanding ( p.5|) in powers 
of the opacity n{^)a{y) p|, The second one, which we will follow in this work^, corresponds to 



(2.8) 



where g(^) is the transport coefficient which characterizes the medium-induced transverse mo- 
mentum squared (g^) transferred to the projectile per unit path length A. Provided the collision 
geometry is known, all medium properties are encoded in the transport coefficient. In the approxi- 
mation ( |2.8[ ), the path integral in ( |2.5D is equivalent to that of a harmonic oscillator and its solution 
is known. The resulting expressions can be found in pT| . 



^Detailed comparisons of both approximations in the massless and massive cases can be found in and ||26 
respectively. 
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For a static medium^, the spectrum of medium-induced radiated gluons is a function of two 
dimensionless variables uj jujc and defined by 



1 k2 

- r 2 2 ^1. 

= , K = — 

2 qL 



(2.9) 



The shape of this function is shown in Fig. |I] for different values of the energy of the radiated gluon. 
Direct comparison with Eq. ( |2.7| ) defines the medium-induced part of the splitting probability 



a, azat 



(2.10) 



where uj = {1 — z)E and = z{l — z)t are taken in the medium-induced gluon radiation spectrum 
- parton masses have been neglected. 



codr^^/[dcodKT 




C0/C0c=10"'' 
03/(0^=10"'' 

C0/C0(;=1 

C0/C0c=10 



K' 



2 10" 



Figure 1: Single-gluon inclusive spectra versus for = 1 = Cji and different values oiuj/tOc- 



At this level, g ^ gg and q ^ qg splitting functions are different by simply the Casimir 
color factors. Let us now discuss finite-x corrections to the leading behavior described in the 
expressions above. The helicity averaged contribution for the case when the radiation vertices (in 
amplitude and complex conjugate amplitude) take place outside the medium recovers the usual 
vacuum splitting functions^ and the full lowest order form will be used for them. In the medium 
case, the helicity averaged contributions to the radiation vertex factorize as a multiplicative factor 



^The case of a medium which expands and gets dilute can be mapped onto a static one with a redefinition of the 
transport coefficient |2S 



29, BO 



■^The reason for this is that in the cifconal approximation, the only effect of the medium on the parent parton 
before the splitting is a color rotation which cancels in the medium averages ^, 0] - see also jsT) for a recent 
discussion. 
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of the form (1 — z)P^^'^{z) / {2Cr) = (1 + z'^)/2 for the case where the parent parton is a quark 
0' S 0' il ^"^^ similarly for the gluon. So, in summary, to take into account the finite-x correction 
to the medium splitting function, we multiply the right-hand side of ( p.lO|) by (1 + z^)/2 if the 
parent parton is a quark and by z if the parent parton is a gluon. In the second case we impose a 
symmetrization around z = 1/2. We have also checked that making the substitution uo — >■ z{l — z)E 
in ( 2.101) 10 produces only small corrections. The impact of these large x modifications is small as 



we show numerically at the level of the Sudakov form factors in the next Subsection. 

For the coupling constant as(Q^), we will use a lowest-order running with 3 quark flavors, 
with Aqcd = 0.236 GeV which gives as{mz) = 0.1172. In the infrared, we freeze the coupling at 
Q"^ = 4 GeV^. The scale is taken to be = z{l — z)t which corresponds to the squared transverse 
momentum of the emission^. This choice corresponds to the implementation of angular ordering 



in the emissions, see e.g. We have checked that all qualitative conclusions we extract with 

this implementation hold when we implement a different scale for the running of the coupling (e.g. 
= t) or when we freeze the coupling to a larger value (e.g. at = 2 GeV^). 

2.2 Sudakov form factors and medium-modified fragmentation functions 

The Sudakov form factor, including summation of different parton species and the relevant kine- 



matical limits reads 22 



Aiit) = exp 



(2.11) 



with Zminit) = to/t. Here, to is some scale below which non-perturbative effects (i.e. hadronization) 
become important. We take to = 1 GeV^. Note that the evolution only makes sense from t 
up to 2to, to prevent Zmin{t) > 1 — Zminif)- At this level, we consider to be independent of 
the parent or daughters being gluons or any quark flavors. In ( p.ll| ) the summation over flavors 



(j = (?, M, M, (i, . . . ) is made explicit to indicate that all possibilities q qg, g gg and g ^ qq 
are taken into account in the evolution. For the case g — > qq, the vacuum splitting function is 
always taken, so it acquires no medium modification as it is subleading in energy; numerically, this 
contribution is anyway very small^. 

Using the medium- modified splitting functions defined by ( |2.4| ), we compute the corresponding 
Sudakov form factors for quarks and gluons. In Fig. ^ we show an example for several values of 
medium parameters q and L, and for parton energies E = 10 and 100 GeV. The influence of the 
medium, and the difference for quarks and gluons, are clearly visible (and the relation of the splitting 
functions at small x given by the Casimir factors almost exactly fulfilled). We also demonstrate the 
small effect produced by the large-x corrections that we introduce. 



^The medium-induced gluon spectrum has been deduced for a fixed coupling constant. Nevertheless, arguments 
have been given |^ that the transverse momentum of the emitted gluon provides the right scale in this problem. 

^Notice that when the evolution takes place in the medium, there is a chance that the parent parton fuses with a 
parton from the medium, gg g, qg ^ q, qq — > g, and conversions of a gluon jet into a quark one or viceversa may 
take place. This possibility has been taken into account in e.g. |3^. The effect of these corrections, not included 
in the present work, is not large and restricted to the small- and moderate-pr part of the spectrum. 
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Figure 2: Sudakov forms factors for quarks (upper curves and symbols in each plot) and gluons (lower 
curves and symbols in each plot), for parton energies S = 10 (plots on the left) and 100 (plots on the right) 
GeV, and for different medium parameters and large-x corrections to the medium splitting functions. In 



the legends on the plots, 'large x' refers to the (default) two first corrections explained in Subsection 2J. 
and 'Zakharov corrections' to the third correction explained in that Subsection. 



In order to compute the medium-modified fragmentation functions, the only missing ingredient 



is the initial condition for the DGLAP evolution ( |2.1| ). This initial condition is a non-perturbative 



quantity describing the formation of the final hadron. Although a medium-modification of this 
quantity could be possible, we will simply assume that this initial condition is unchanged and equal 
to the vacuum one 

D"^^'^(x,to)=I^""^(x,to), (2.12) 
and we take D™'^(x, to) from ||3^. The motivation for this ansatz is the following: in hadronic 



collisions, particles produced at high enough transverse momentum hadronize outside the medium. 
Eq. (p. 12 ) assumes that this non-perturbative hadronization is not modified by the medium, whose 
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effect is only to modify tlie perturbative associated radiation^. All present radiative energy loss 
formalisms rely on this assumption (see the Appendix ^ for the equivalence with these formalisms). 

The medium- modified evolution given by Eqs. (|2.2|) , ( p.3|) and ( |2.4| ) and the initial condition 
( p.l2| ) complete the formulation of our approach. 



3. Results 

We solve Eq. (|2.2|) numerically in order to obtain our medium-modified fragmentation functions 
(MMFF). We use a brute-force method with a 4th-order adaptative Runge-Kutta for the evolution 
in virtuality plus a Gaussian quadrature for the integration in momentum fraction (and virtuality 
for the Sudakov form factors). For the initial conditions we use the Kniehl-Kramer-Potter (KKP) 
p3| parameterizations aX = 2 GeV^ for pions and we consider only three flavors u, d, s^. We 
have checked that our simplified evolution reproduces the KKP results better than 40 % in the 
region of validity of the KKP parametrizations (0.1 < z < 0.9) up to the highest virtualities we 
have considered. 

In Figs. 1^ and ^ we show our results for the medium-evolved fragmentation functions onto pions 
for two different parton energies (10 GeV, relevant for RHIC and 100 GeV, relevant for the LHC), for 
different medium parameters and different parton types. The main effect found is a suppression of 
the fragmentation functions at large z values and a corresponding enhancement at small z. These 
two features grow with increasing scale (the higher the scale, the higher the length of evolution 
thus the greater the medium modifications). They also increase with the transport coefficient and, 
as shown in Fig. with the medium length. These medium effects decrease however with the 
parent parton energy. The qualitative reason for this is that the energy loss becomes more and 
more energy-independent with increasing energy and thus the fractional energy loss decreases with 
increasing energy. As expected, the medium modifications are greater for gluons than for quarks. 

In order to check the feasibility of our approach and the compatibility of the extracted medium 
parameters with those obtained within other approaches, we compute the medium-modified spectra 
of neutral pions at RHIC. We standardly convolute our medium-modified fragmentation functions 
with the nuclear parton densities (we use the EKS98 parametrization from Ref. ||3^) and the 
hard scattering elements to obtain the medium-modified particle spectra following the factorization 
formula at LO: 



a 



fAixi, Q^)fB{x2, Q^) ® a(xi, X2, g') ® A^h(^, Q^). (3.1) 



As it is common in the phenomenology of jet quenching in heavy ion collisions, we show our results 
for the nuclear modification factor defined as 



^For processes with different kinematic conditions |Q this assumption could not hold, but the negligible effects 
seen in dAu data at RHIC indicate that this is a reasonable assumption for particle production at high-pt in 
nuclear collisions. 

^We have checked that the use of a different parametrization |^ for the initial conditions does not change the 
conclusions of our study. 
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Figure 3: Medium-modified fragmentation functions for gluons (left), u (or d) quarks (middle) and s 
quarks (right) as a function of the energy fraction carried by the hadron, z. The parent parton energy is 
10 GeV and the medium length is L = 2 fm. For each of the three different scales (Q^ = 2, 19 and 400 
GeV^ in black, green and red respectively) we plot the fragmentation functions at three different medium 
densities: vacuum (solid), q = \ GeV^/fm (dashed) and g = 10 GeV^/fm (dotted). 



where the numerator correspond to the cross section ( |3.1| ) including the nuclear corrections to the 
PDFs and our MMFF and the denominator is computed with the proton PDFs and the vacuum 
fragmentation functions. 

In order to check the validity of our formalism we first take a fixed length of L = 6 fm ~ Rku 
in the calculation of the MMFF. This simple geometry favors a value of g ~ 1 GeV^/fm - see Fig. 
^- in agreement with the findings in [rT|, where also a fixed length was taken. As it is known, the 
actual geometry used in the calculation affects the extracted value of the transport coefficient when 
the experimental data is fitted. It is not the goal of this paper to repeat fits of experimental data 
which have been extensively discussed in the literature in the past years. Here, we simply repeat the 
calculation with a geometry which takes into account the production point of the high-p^ parton 
inside the extended medium. The procedure, performed for a cylindrical (constant profile) and for 
a spherical (constant density) nucleus, is as follows: i) we first sample a production point and a 
emission angle inside the medium by Monte Carlo; ii) the traversed length is then calculated and 
used to compute the cross section (|3.1|) with medium effects; iii) the points i) and ii) are repeated 
to obtain an average cross section in nuclear collisions; iv) the ratio ( p.2| ) is done. The obtained 
result is plotted in Fig. ^ where the favored value of q is larger by roughly one order of magnitude. 
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Figure 4: The same as Fig. |3| but for a parent parton energy of Ejet = 100 GeV and scales = 2, 137 
and 40000 GeV^ in black, green and red respectively. 



Similar results were obtained in [3S 



As a last comment, we have set the fragmentation scales equal to the fragmenting parton 
transverse momentum. Other common choice is to set the factorization scale to the final hadron 
transverse momentum, which is smaller by a factor z. In Fig. ^we show that both give similar results 
for the proton-proton case once the K-factors are adjusted to describe the data. However, the second 
choice would reduce the suppression, as it is evident from the Figs. || and |. Although the choice 
of scale is something arbitrary, this fact demonstrates the role of virtuality effects in jet quenching 
calculations. We find our implementation more natural in the picture of a parton branching process 
where the total amount of radiation is determined by the initial maximum virtuality. The latter is 
dictated by the perturbative hard cross section. 



4. Comparison with the approach based on quenching weights 



Medium modified fragmentation functions are usually computed by an energy shift of the corre- 
sponding vacuum fragmentation functions P0|: 



de 



Pie) D' 



X 



(4.1) 
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Figure 5: Left: Fragmentation function for gluons onto pions at Ejet = 100 GeV. Our results are plotted 
at = Ej^^, for three different mediums: vacuum (black), q = 10 GeV^/fm (green) and q = 50 GeV^/fm 
(red) and for two different medium lengths: 2 fm (solid) and 6 fm (dashed). Right: Medium to vacuum 
ratio of the gluon fragmentation functions for the same values as in the plot on the left. 
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Figure 6: Left: Nuclear modification factor Raa computed with the obtained medium- modified fragmen- 
tation functions for a fixed in-medium path-lenght of L = 6 fm. Right: Same but computed with more 
realistic geometries leading to a distribution of path-lenghts over which the suppression is averaged. In 
both cases, the experimental data is taken from |3q|. 



This model assumes that a highly energetic parton losses a fractional amount of energy e while 
traveling through the medium and fragments with un-modified (vacuum) fragmentation functions 
once it is outside. Any modification of the virtuality dependence of the fragmentation is neglected 
and the probability distribution for the energy losses - quenching weights - has a discrete and a 
continuous part, 



(4.2) 



given by 



Po = exp 
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(4.3) 



(4.4) 



In Appendix ^ we show that this procedure corresponds to the hmiting case of our formahsm for 
very large virtuahties and considering only soft emissions. Here we show in Fig. [7| the comparison 
of the results of our medium- modified evolution with those using the quenching weights []1T|. It 



can be seen that with increasing length of the evolution, the results of both approaches get closer 
and closer - specially for large values of z most relevant for phenomenological applications -, in 
agreement with the discussion in Appendix 0. 



Dashed - Quenching weights 
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Figure 7: Left: Fragmentation function for gluons onto pions computed with our medium-modified 
evolution (solid lines) and through the standard convolution with the quenching weights (dashed lines), 
for Ejet = 40 GeV, q = 1 GeVVfm, L = 6 fm, and different = 2, 300 and 1600 GeV^. Right: Ratio of 
the fragmentation function for gluons onto pions in a medium with the same characteristic as for the plot 
on the left, over the fragmentation function in the vacuum, for the same values of Q^. 



5. Conclusions 

We have presented a new implementation of medium effects in shower evolution which takes into 
account vacuum and medium-induced splittings on the same footing. We introduce a medium- 
modified splitting function into a DGLAP-like evolution. By doing so, our formalism improves 
previous ones, based on the Poisson approximation for the multiple gluon emission, by a consistent 
implementation of kinematic constrains on every individual splitting as well as on secondary emis- 
sions, which are included automatically . The explicit dependence of the virtuality in the evolution 
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makes it possible to study its effects for tlie first time in tlie case of medium-induced gluon radi- 
ation. As expected, we obtain a softening of tlie fragmentation function. This softening increases 
with increasing virtuality, and with increasing medium length and density. In our implementation 
any effect on the initial non-perturbative condition is neglected since we do not expect it to change 
if large enough virtualities and parton energies are studied. The limited momentum range of ap- 
plicability of the DGLAP fits |33] makes not possible to predict with reliability the crossing point 
between parton suppression at large-z and parton enhancement at small-z, which also depends on 
the initial conditions chosen for the evolution |]33|, |35|1. This limitation will be circumvented in future 
works. 

One of the interests of our formalism is the simplicity to be included into a Monte Carlo code 
by implementing a medium-induced splitting probability AP{z, t) into the parton shower evolution 
routines. Work along this direction is in progress. In particular, it will allow to compute observables 
beyond single inclusive production, like two-and three-particle correlations, and compare the results 
with other approaches It will also allow to study the possible interplay between the finite 

formation time of the partons radiated in subsequent emissions and the length of the medium. 



which has been ignored here as it is ignored in existing formalisms of energy loss [11, 12, 13, 14 



Finally, we have provided an independent check on the quality of the Poisson approximation 
for the calculation of inclusive particle production at high px- We have presented analytical and 
numerical checks of this approximation, and shown the agreement between both approaches for large 
values of the fraction of momentum. These are precisely the relevant values for the calculation of 
the inclusive particle suppression in heavy-ion collisions. This provides a connection with all the 
previous phenomenology and further supports the usual assumptions. 
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A. From the Sudakov form factors to the quenching weights 

In this Appendix we show how the quenching weights []TT], |T^ appear when considering the DGLAP 
evolution of the fragmentation functions. Ignoring for simplicity all parton and particles labels. 
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DGLAP evolution can be written as 



D{x, t) = A{t)D{x, to) + A{t) / V (? ^0 • ^^-^^ 

The iterative solution to this equation reads 

D{x,t) = Ait)Dix,to) + A{t)f] r^... ['"'^^ 1^1 I ^ 

~^Jto '^l J to ^2 J to tn J Zi J Z2 J Zn 

X Piz,)P{z,) ■ ■ ■ P{z^)D ( , to 

\Z1Z2 ■ ■ ■ Zn 

= Ait)Dix, to) + A(t) / f r ^ r ^ . . . '^fifdz. p(.. 



to ^1 Jto ^-2 Jto .^-^ 



X 5 {ziZ2 ■■■Zn-[l-t])D toj . (A.2) 

Now we consider = 1 — <C 1, j = 1, 2, . . . , n i.e. the successive emissions to take a small 
fraction of the parent energy-momentum, which results in 

5 {Z^Z2 ■■■Zn-[l-e])c^5i^ " ^''') ^^'^^ 



and 



D[x,t) A{t)D{xM) + m I ^ [ dz^P{z,) 

J ^ ^ n=l i=l -Jio J 

Equation ( |A.4|) holds both for emissions in the vacuum (with A™'^(t), P™'^(z) and D^'^^{x,t)) 
and in the medium (with A™'^'^(t), AP{z) and D^'^'^{x, t)). Taking into account that in our approach, 
due to the property (|2.6|), we have 



P{z) = P™=(z) + AP{z), A{t) = A™^(t)A™<='^(t), (A.5) 

we get 



D{x,t) ~ A^^'i(t)D™'=(x,t) + A^<='^(t) [ -^J^l-f[ f ^ f dziAP{zi) 

i 1 ^ n=l i=l J 

(A.6) 
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where the sum and product run now over medium-induced emissions and we have used the fact that 



\ vac mcd / 



X 



(A.7) 



mod 




X 



e')(l 



with d-e'-e" 



(1 -e')(l -e") for e',e" < 1. 
Finally, we ignore virtuality (as done in all previous approaches to radiative energy loss except 
in [^) and set in the integrals over tj and Zi the kinematical limits for gluon emission. Then we can 
identify, using (|2.9|), (|2.7| ) and (p.lO|) , the discrete and continuous parts of the quenching weights 
lTT|,IT2i 



Po = exp 



dcu / (ik_L 



dl 



med 



dudk 



oo n „ „ 
n=l i=l 



dl 



mod 



duidh^i 




(A.8) 
(A.9) 



to get the standard expression for the fragmentation functions in a medium through the quenching 
weights pOj : 

D{x,t)c:,poD-^'^{x,t) + J ^p(e)D™^ (r^'*^) ■ ^^-'^^^ 

Summarizing, DGLAP evolution of fragmentation functions in a medium corresponds to the 
usual picture of radiative energy loss through quenching weights when virtualities are ignored, and 
only soft emissions are considered. This latter consideration lies at the root of all existing formalisms 
for radiative energy loss. In this correspondence, the discrete part of the quenching weights maps 
onto the Sudakov form factor, ( |A.8| ). 
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